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Abstract 

The quantum two-center MICZ-Kepler system is considered in the limit when one of the interaction centers is 
situated at infinity, which leads to homogeneous electric and magnetic fields appearing in the system. The emerging 
system admits separation of variables in the Schrodinger equation and is integrable at the classical level. The first 
order corrections to the unperturbed spectrum of the ordinary MICZ-Kepler system are calculated. Particularly, 
the linear Zeeman-effect and effects of MICZ-terms are analyzed. The possible realizations of the system in some 
quantum dots are considered. 
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1 Introduction 

The MICZ-Kepler system describes an electrically charged scalar particle moving in the field of static Dirac dyon(s), 
i.e. a particle carrying both electrical and magnetic charges. This system was originally proposed by Zwanziger 
P] and Mcintosh and Cisneros independently [2 . In their papers the corresponding system was obtained from the 
ordinary Coulomb problem by replacing the interaction center by a dyon which has, beside the scalar potential, the 
vector-potential as well, and by introducing into the Hamiltonian the additional centrifugal-like term, s 2 /2/ir 2 , where 
s = eg is the monopole number (we use c = h = 1 units), e is the electric charge of the probe particle, g is the 
magnetic charge of the static dyon and [i is the mass of the probe particle. According to the Dirac quantization 
rule the monopole number can take integer or half-integer values: s =, 0, ±1/2, ±1... The corresponding system is 
described by the following Hamiltonian: 
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Hmicz = — (p - eA g ) 2 + ^—^ + — , (1) 



rotA g = (2) 



Due to this additional potential term the emerging system inherits all distinguished features of the underlying Coulomb 
system. For instance, the MICZ-Kepler system, beside the general rotational symmetry generated by the algebra of 
angular momentum components (so(3)), has also a higher so(4)-hidden symmetry which is connected to the conser- 
vation of the Laplace-Runge-Lentz vector. The expressions for these conserved quantities are very similar to those of 
the Coulomb problem 

r 1 r 

J = r x (p — eA) + s-, I = — (p — eA) x J — eq-. (3) 

r fj, r 

As one can see, only the additional spin-like term appears in the expression for the orbital momentum. The hidden 
symmetry of the MICZ-Kepler system leads to the possibility of separation of variables in several coordinate systems, e. 
g. in spherical and parabolic ones [3]. The origin of the additional potential term (MICZ-term) could be interpreted as 
the interaction energy between the spherically symmetric magnetic field of Dirac monopole and the magnetic moment 
of the probe particle, M. = <^;J, where J is given by Eq. ||3J). However, additional part of the angular momentum 
vector proportional to the monopole number s should be rather assigned to the monopole electromagnetic field by itself 
but not to the particle. Thus, such a clear and natural, at the first glance, explanation of the origin of MICZ-term is 
not quite correct. On the other hand, this term can be obtained both on classical and quantum level when one reduce 
the four-dimensional isotropic oscillator with the aid of so-called Kustaanheimo-Stiefel transformations [4] . 

The shape of the classical trajectories of the MICZ-Kepler system coincides with these of the ordinary Coulomb 
one, however in contrast to the latter the orbital plane are non orthogonal to the angular momentum vector. The 
angle between J and orbital plane normal satisfy the relation: 

cosx =w (4) 



1 



which immediately followed from Eq. |3]). The properties of the solutions of the quantum mechanical MICZ-Kepler 
problem undergo minor modifications with respect to those of the underlying Coulomb one. Namely, it has the same 
spectrum, given in terms of the same quantum numbers. The only difference is in the range of the possible values of 
angular momentum. The presence of the monopole with monopole number s shifts upward the allowed values of orbital 
and magnetic quantum numbers j and m: j = |s|,|s|+l,|s|+2,... There is also more general statement concerning 
the incorporation of Dirac monopoles in spherically symmetric mechanical systems. Any spherically symmetric system 
(without monopoles) defined on the conformally flat space with metric G(r) will preserve its spectral properties except 
the possible range of angular momentum values mentioned above, if adding the monopole potential will accompany 
with the following modification of potential term [U [6l [JJ [8] : 

U(r) - U(r) + (5) 

Another difference between Coulomb systems and their MICZ-counterparts manifests itself in the behavior in external 
fields. So, one can observe the modification of the selection rules in the dipole transitions J9[[l0]. The investigations 
of the generalizations of the MICZ-Kepler system to the higher dimensions and/or curved spaces as well revealed 
their close similarity to the underlying Coulomb systems [7J QTJ [TJ1 [T3] . Multi-center generalization of MICZ-Kepler 
systems has been found recently in Ref. pj Q3] . The Hamiltonian which describes the charged particle moving in the 
electromagnetic field of n Dirac dyons fixed in the Euclidian space points with radius- vectors a^, i — 1, n reads: 

1 / ™ \ 2 n 2 / " \ 2 

*-£(p-«|>.C'-.>J +-E F ^ r +s:feiF^i) • «**«m-£ < 6 > 

where qi(gt) are the electric(magnetic) charges of the the i-th dyon. The last term in the expression generalizes the 
MICZ-term for the multi-center situation. In full analogy with the corresponding pure Coulomb system two-center 
case is integrable on the classical level [Jj [14] , i .e. it allows separation of variables in the Hamilton-Jacobi equation 
in elliptic coordinates. The limiting case of the two-center MICZ-Kepler system when one of the dyons is situated 
at infinity which results in the homogeneous electric and magnetic fields, is also integrable on classical level. The 
MICZ-term corresponding to this dyon results in its turn in the quadratic potential. Corresponding Hamilton-Jacobi 
equation become separable in parabolic coordinates. Because of presence of homogeneous electric and magnetic fields 
imposed to the one-center MICZ-Kepler (charge-dyon) system one can call it MICZ-Kepler-Zeeman-Stark system. 
In this paper we consider some quantum mechanical issues of the MICZ-Kepler-Zeeman-Stark system. In addition to 
the pure stark effect in charge-dyon system considered earlier in Refs.[T"5l 110] . we present perturbative corrections to 
the spectrum for the charge-dyon system corresponding to the homogeneous magnetic field(linear Zeeman effect). 

It is also noteworthy that multi-center MICZ-Kepler system allows Af = 4 supersymmetric extension when mag- 
netic and electric charges of all dyons satisfy the trivial Dirac-Schwinger-Zwanziger condition 14J 

9iQj - 9jH = 0- (7) 

The corresponding Af — 4 supersymmetric mechanical systems were constructed in Refs. [HI [UJ and with a little 
different approach leading to the MICZ-Kepler system on the three-dimensional sphere in Ref. [T5]F1 

The paper is organized as follows. In the second section we derive the Schrodinger equation for the MICZ-Kepler- 
Zeeman-Stark system in parabolic coordinates. The solution of the unperturbed problem is presented and calculated 
the first and second order corrections in B and £ The average dipole and magnetic moments acquired the system are 
also calculated. In the third section the related quantum dot models are discussed. Some final remarks are presented 
in the conclusion. The appendix contains some technical points. 

2 Two-center MICZ-Kepler system and quantum MICZ-Kepler— Stark— 
Zeeman System 

The Hamiltonian of the two center MICZ-Kepler system on Euclidean space has the form [T?] 

n =\ (p - eA - (r - ri) - eA - (r - r2))2 + \ (f^f + ra) 2 + e (k^i + ra) ' (8) 



x For a unified superfield formulation of N=4 off-shell supermultiplets in one space-time dimension see | 19| . For earlier work see 1201 . 
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where 

A g (r)=g— 9 

r(r — nr) 

is the vector-potential of the Dirac dyon with magnetic charge g, n is the unit vector pointing to the singularity line. 
Choosing an appropriate gauge for the monopole vector-potential ([9]), 

A r =A e = 0, A v = g cos 9, (10) 

and supposing that the dyons be fixed on the z-axis at points (0,0, a) and (0,0,— a) respectively, one arrives at the 
following expression for the Hamiltonian (jHJl in spherical coordinates: 

H4^4+ ^" lC f 1 n 2C0Sg2) V !f- + -) 2 + - + - > (11) 
2 r r 2 r 2 sin 2 9 2 \n r 2 ) r x r 2 V ; 



where r\ ; 2 — \J x 2 + y 2 + (z ± a) 2 . As it was shown in Rcf. JT5] this system on the classical level admits separation of 
variables in elliptic coordinates, what leads to the integrable system of Hamilton-Jacobi equations. The two-center 
MICZ-Kepler system has an important limiting case, when one of the interaction centers is situated at infinity. The 
field of such a dyon results in the homogeneous electric and magnetic fields being parallel to each other. If one takes 
for the vector potential of the homogeneous magnetic field the following gauge: 

A r =A 9 = 0, A v = X -Br 2 sin 2 9, (12) 

then the corresponding Hamiltonian in spherical coordinates reads 

f 

r 2 r 2 sin 2 9 J 2 \r / r 

where B and £ stand for the modulus of the magnetic and electric field respectively which are pointed in the z 
direction. The classical case is characterized by the separation of variables in the Hamilton-Jacobi equation in 
parabolic coordinates given by the following relations: 

y 

£ = r + z, r/~r~z, </? = arctan — . (14) 

x 



n, 1/ 2 pj (p v -scos9-±eBr 2 sm 2 9) \ 1 fs n \ 2 eg 

H = o [Pr + ~2 H , . 2n + o ~ + eBz H e£ Z > ( 13 ) 



The same feature holds in the quantum case for the Schrodinger equation for the Hamiltonian (| 13|) as well. The 
quantum Hamiltonian resulting from Eq. (|13p reads 

1a / scos9 \ d s 2 , , , 9 > eq 

H = --A+i[ - . 2q +lo b ) — + , + -2- (p 2 + Az 2 ) + Zslub cos9 + -1 - e£z, 15 
2 \r 2 sm z 9 ) dtp 2r 2 snr0 2 x ' r 

where ub = is the cyclotron frequency in our units (e = K = fj,= 1) and p 2 = r 2 sin 2 9 — x 2 + y 2 . Being rewritten 
in the parabolic coordinates given by Eqs. (fTl)) , the corresponding Schrodinger equation takes the form 

d A <9*\ d ( dm\\ 1 f d 2 ,\ , „. / s (\ 1\ \ 9* 



i + V + d~V W ) + Tr, \W 2 - 8 ) * + 21 [i^J U ~V ) " ) ^ + (16) 



2 £ - LZ. _ 6flCl , * - _2- + e f ^ — * = 0. 

Hereafter in order to deal with discrete spectrum we suppose the electric charge e of the probe particle to be negative. 
Seeking for the eigenfunctions in the form 



where m is the magnetic quantum number, and dividing the equation on 4 /l^ 2 , one obtains 



* = M0Mv)e imv , (17) 

1/1/2 

t+Tl 
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where 

u ® = ~^- !! T? + T? + w -t' (19) 
v(v) = - s ^ ) - l fv 3 - e 4v 2 + w +v 

and the following notations are introduced s± — m ± s, W± — 1/4(215 + 2mu>B ± 6so->b). Separating the variables £ 
and 77 one obtains the system of equation for the functions /1 and / 2 

U (d i) + (^'-i-^ 3+ ^ e+w - ( ) h =°- (20) 



drj \ drj J \ 4rj 4 4 

where the separation constants /?i and /3 2 satisfy the condition 

Pi +02=1 (21) 
2.1 One-center MICZ— Kepler system in parabolic coordinates 

Note that the case of the one-center MICZ-Kepler (B = 0, £ = 0) system completely coincides with the quantum- 
mechanical problem of the hydrogen atom. It leads to the same spectrum. The only difference consists in the 
allowed values of magnetic quantum number, which are just shifted by s. Now one can regard the terms corresponding 
to the homogeneous electric and magnetic fields as perturbations with respect to the integrable case of the one center 
MICZ-Kepler Hamiltonian which leads to the following equations[3]: 

i(^) + (l^- S -k + \e\ q Pi)f^0, (22) 



d£ V d£ J V 2 4 £ 

J;("f) + (^-f + '*A->=° 

which are almost identical to those for the hydrogen atom in parabolic coordinates (see for example [22]). We consider 
only the discrete spectrum of the unperturbed problem (E < 0). In this case it is convenient to pass to the new 
variables n = l/\/—2E 1 pi = £/n and p 2 — rj/n. Let us write down the equation for the function /1 in these 
coordinates 

d2 ^ m 1 dfl x ( 1 -l 1 ( jl 1 n I xi0l s + \t n 
dpi pi dpi V 4 Pi V 2 J Apf J 

where we introduce the notations and put for simplicity \e\q = 1 (natural units) 

m = npi - 1 (|s+| + 1) , n 2 =n/3 2 -i(|s_| + l). (24) 

The corresponding solution for / 2 can be obtained by replacing Pi by P% and s+ by s_. Taking into account the long- 
and short-distance asymptotic behavior of the solutions one can seek the function fi in the following form: [22j 

h(p 1 ) = e-^ 2 p[ s+W2 F 1 ( Pl ), (25) 
which leads to the equation for the confluent hypergeometric function for F± (pi ) 



d 2 F^ dFi 
Pi^rr + (\8+\ + l-p 1 )=±+n 1 F 1 =0. (26) 
dpf dpi 



Thus, the functions F± and F2 which satisfy the finiteness condition read 

F 1 (p 1 ) = F(-n 1 ,|s + | + l;p 1 ), F 2 {pi) = F (-n 2 , |«_| + l;p 2 ) , (27) 



2 Related supersymmetric systems with Dirac monopoles as well as the details of the Hamiltonian reduction technique have been 
considered in 1211. 
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where n\ and n 2 must be non-negative integers. So, from Eq. (124|) one obtains that the principal quantum number 
for the one-center MICZ-Kepler system is 



m + s + m — s , 
n = m+7i 2 + J L_J L + l ) (28) 

Thus, for the unperturbed part of the Hamiltonian (|13p . each stationary state of the discrete spectrum in parabolic 
coordinates is characterized by three integer quantum numbers, i.e. the parabolic quantum numbers ri\ and n 2 and 
the magnetic quantum number m. The normalized eigenfunctions are |22j 

Vnuw&w) = -^7=fn 1 ,r n +.(Z/n)f na ,m-.(S/n)j mv , (29) 

where 



fk * {p) = p\/™^V pl ^ (-*. \p\ + 1; p) ■ (30) 

2.2 The perturbative first order corrections and linear Zeeman splitting 

One can regard all terms, corresponding to magnetic and electric held as small perturbations. Thus, we have the 
unperturbed Hamiltonian of the one-center MICZ-Kepler system with exactly known stationary state eigenfunctions 
(j2"9"|) and the perturbation given by the operator 

W . W . + + 4,*) _ If (f - „ - + *f*Lg. (3!) 

There is also a constant shift in the energy level, caused by the magnetic field and proportional to the magnetic quantum 
number m, i.e. —mu>B- So, the first order correction to the stationary unperturbed hydrogen-like spectrum 3 

E nin 2 m = } j ~7~ T2> (32) 

2 I ni + Tl2 + 1 + J L 2 J 1 J 

reads 

(ni,n 2 ,m|yV|ni,n 2 ,m), (33) 

Thus 

E^Lm = Ei nL m - mw B + AEW n2m , (34) 

where 



, m \e\£^ 3s|e|R ne 2 B 2 , 
AE 1 ' m = — X" IT H 35) 



and it stand for the integrals of hypergeometric functions presented in Appendix. Here we restrict ourselves with 
the first-order corrections in £ and B (we will suppose £ and B to be of the same order). Thus, 

E £Lm = -^2 + |e|£(nn_ - ™) - ^^(^ - m) + 0(S 2 ), (36) 

where 

|m + s| — I to — s| 



+ ni — n 2 . (37) 

Putting B = one obtains the well known result of purely Stark-effect in the charge-dyon system (TU1 [TS] . The third 
term in the expression corresponds to the linear Zeeman-effect. As usual, in this approximation magnetic field removes 
the degeneracy with respect to the magnetic quantum number m and induce the magnetic moment in the system. The 
average magnitude of induced magnetic moment could be obtained by taking the derivative of Eq. (|36p with respect 
to B. However, the term 3s ^ B cos 9 in the operator (f3"Tj) actually comprise of two parts of different origin. The first 
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is 2 cosO which is originated from the magnetic field in the kinetic term. The second part is s\e\B cos 9 coming 
from the additional potential term (MICZ-term). Indeed, the MICZ-term in this system just has the coefficient which 
is coincide with magnetic field by magnitude but physically is not identical with it. Thus, in order to calculate the 
induced magnetic moment of the charge-dyon system in the external magnetic field one must subtract from En^ n2m 
the impact of the MICZ-term prior to taking the derivative: 



,S7l- 



M z = ^ = -|e|(— - m) (38) 

( \m + s\ — \m - 

= Ms 



\Enin 2 m ~ s\e\B(ni, n 2 , m\ cos 9\ni , n 2l m)^j i 
dB 

|m + s\ — \m — s\ + 2(ni — n 2 ) 
\m + s\ + \m — s\ + 2(ni + n 2 + 1) 

where hb = \\&\ is the Bohr magneton in our units. Taking into account above mentioned properties of the linear in 
B corrections one can write down the spectrum in the following form: 

E$ n2m = 4°Lm -d z S-M z B + . (39) 

The last term originated from the MICZ-term and 

— d/\Enin 2 m 3. . . TOS . . , 

d£ = ~2~' e ' f ( nn -~~3")' ( 40 ) 

is the average dipole momentum in linear approximation, calculated in Ref. [15] For the ground state of unperturbed 
system (charge-dyon system) which is characterized by n\ — n 2 — 0, |m| < \s\ n = \s\ + 1, n_ = msgn(s), 
m = — \s\,— \s\ + 1,..., \s\ — l,\s\ one obtains 

m 

M z = -M BTTR . (41) 

Thus, the linear Zeeman-effect in the charge-dyon system removes the (2|s| + l)-fold degeneracy by m of the ground 
state. The induced magnetic momentum for positive values of m is always negative and goes to zero at \s\ — *■ oo. At 
s = one get the expected value A4 Z = — /Xf,m corresponding to the ordinary Zeeman splitting in the hydrogen atom. 
Hence, for the ground state of the system under consideration one can write 

E «' - -idw + m|e|sg " s ( w + 1) £ + ^TWl 8 + |e| irw a m 

As one can see, the sign of the linear Zeeman effect depends only on the sign of the magnetic quantum number m, 
whereas the linear Stark effect depends on the relative sign of m and s. If we suggest that coefficient in the oscillatory 

2 

potential -Yip 2 + 4z 2 ) in Eq. ([3T]) is independent of B and is of the same order with £ and B then according to the 
formula ([33]) we will get the corresponding first order correction to the ground state energy in the following form (see 
Appendix) 

AE^ = y(l + N)((»i 2 + s 2 )(\s\ + 6) + |.s|(2m 2 + 11) + 6). (43) 

3 Additional potential term and related quantum dot models 

It is obviously seen that one can add to the Hamiltonian (|13[) an additional potential term of the form 

without breaking its classical integrability and separation of variables in the quantum case in parabolic coordinates. 
Such kind of potential can originate in the model of a cylindrical quantum dot, where the influence of the dot boundary 
is described by the confining potential of the parabolic type (oscillator potential, see also [23]) with special rate of 
radial and transverse frequencies. Moreover, if we put s = the emerging system could be identified with the model 
of a charged particle moving in the axially symmetric quantum dot in the field of Coulomb center and homogeneous 
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electric and magnetic field pointing along the dot symmetry axisH The corresponding confining potential takes the 
form 

f 



U c = - (lu 2 P 2 + Cj 2 z 2 ) , u> = 2^uj 2 +lu 2 b . (45) 

This model is integrable at the classical level [Mj [7] and leads to the separable Schrodinger equation. One can also 
consider a little different context of that issue. Let us suppose we have a particle moving in the field of Coulomb center 
(or the relative motion of two electrically charged particles) in the cylindrical quantum dot with an axially symmetric 
oscillatory confining potential with arbitrary frequencies uj p and uj z . The Schrodinger equation for such a system does 
not admit separation of variables, except for the case uj z — Auj p . Let us assume also that there are homogeneous 
electric and magnetic fields both pointing along the z-axis 

-e£z + ±(Lo 2 pP 2 +La 2 z z 2 ). (46) 




However, if the magnetic field magnitude B takes the value iyw^ — Aio 2 p the variables separate in parabolic coordinate. 

So, at the classical level one obtains an integrable system, which can be a subject for the perturbation theory in the 
quantum case. The quantum dots models with Dirac monopole inside has not only academician interest. Though, 
no elementary particle carrying magnetic charge has been discovered up to now, there are very promising theoretical 
evidences of the magnetic monopoles existence as emergent particles, i.e., as the quasi-particles in some strongly 
correlated many-body systems. It was shown in Ref. [25) that such kind of magnetic monopoles do emerge in so- 
called spin-ice materials, the exotic class of magnets in which local magnetic moments (spins) are residing on the 
sites of pyrochlore lattice (corner shared tetrahedra) and are constrained to point along their local Ising axes. The 
spins interact to each other via nearest neighbor exchange and long-range dipole-dipole interaction. Thus, the model 
discussing in this section can be realized to certain extent in the quantum dots prepared from spin-ice compound, such 
as Ho2Ti2C>7, Dy2Ti207, e. t. c. (for review of spin-ice see [26]). 

In order to apply perturbative results from previous section to the quantum dots model discussing here the pa- 
rameters of confining potentials must be small(more precisely, the corresponding correction AE nin2Tn must be much 
smaller than the distance between corresponding neighbor levels of unperturbed problem) , otherwise only the numer- 
ical calculations are relevant. The corresponding correction to the ground state energy has the same form as Eq. (j4"5|) 
with ujq = tug + lo. In this case, in order for the perturbative calculation to be relevant, this correction must be much 
smaller than the ground state energy, which leads to 

lj « ((1 + |s|) 3 ((m 2 + s 2 )(\s\ + 6) + |s|(2m 2 + 11) + 6))~ 1/2 . (47) 

A similar model of the two-electron quantum dot subjected to the external homogeneous electric field was considered 
numerically in Ref. [27]. 



4 Conclusion 

In this paper we considered the quantum mechanical two-center MICZ-Kepler system in the limit when one of the 
dyons is situated at infinity.The electro-magnetic field of such a dyon results in the homogeneous electric and magnetic 
fields parallel to each other. At the same time, the specific additional potential term (MICZ-term) transforms into 
oscillator potential in the direction of the external fields and potential proportional to cos 9. Thus, the system 
considered in the present paper concerns the "monopolic" generalization of the hydrogen atom (Kepler problem) 
subjected simultaneously to the constant uniform electric and magnetic fields. We separated the variables in the 
corresponding Shrodinger equation in the parabolic coordinates and analyzed the simplest case of the emergent MICZ- 
Kepler-Stark-Zeeman system when external electric and magnetic fields can be regarded as small perturbations. The 
separation of variables is possible in virtue of the additional potential term. As a starting point we considered the 
Hamiltonian of one-center MICZ-Kepler system with well known exact solution, whereas all other terms were treated 
as perturbations. The exact wave functions in parabolic coordinates were used to develop first order perturbative 
calculations. We obtained linear in £ and B corrections to the unperturbed spectrum of the ordinary MICZ-Kepler 
problem. The Stark-effect in the charge-dyon system was calculated within the perturbation theory up to second 
order earlier in series of paper 15, 10J. Here we investigated the first order corrections corresponding to magnetic field 

3 In the noncommutative framework the Coulombic monopole has been considered in 1241 . 
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effect (linear Zeeman effect). We found the average magnetic moment acquired by the charge-dyon system in this case. 
For the ground state the induced magnetic moment is even function of the monopole number s, M. z {— s) = A4 z (s) 
and vanishes at s — > ±00. It is also non-analytic at s = 0. As usually, magnetic field removes the degeneracy of the 
ground state . The impact of the additional potential, so-called MICZ-term, in the first order has much in common 
with that of the magnetic field. The sign of the corresponding correction depends only on the sign of magnetic 
quantum number m. We also analyzed some condensed matter models which can be relevant in the context of the 
quantum mechanical system considered in the papers [14[ 127] . The recent results concerning possibility of formation 
of the magnetic monopoles as quasi-particles in some strongly correlated spin systems, so-called spin-ice materials 
|25] make the attempts of understanding the behavior of the various (quantum-)mechanical systems at the monopole 
background very important and promising for future experiments. Unfortunately, in order to obtain reliable results in 
the quantum dots models with monopoles in the external fields one can not restrict himself with integrable cases and/or 
first order perturbative calculations. One of the reasons for it is the essential role of the confining potentials which, 
generally speaking, can not be regarded as perturbation. However, the second order corrections in the MICZ-Kepler- 
Stark-Zeeman system can lead to new interesting features, for instance, to the interplay between electric and magnetic 
properties of the charge-dyon system. Namely, the correction to the spectrum proportional to £ B which appears in 
the second order perturbative calculations means that average dipole (magnetic) moment depend on magnetic (electric) 
field. Another important issue is to understand spin effect and spin-orbital coupling in the systems with monopoles 
within the quantum dots models. We intend to consider all these questions in the forthcoming papers. 



5 Appendix 

Here we calculate the general integrals which emerge in the perturbation theory calculations for the quantum mechan- 
ical systems in parabolic coordinates 

/>oc />oo 

Z|=/ / (p k 1 ±p k 2 )ft, m+s (pi)fl, m - s (p2)dp 1 dp 2 (48) 
Jo JO 

Let us plug into the integral the expressions for functions f n;P from Eq. (f30|) . Then 

= C ni m+s C n2 m _ s (J nii+ J n2 _ ± J„ 2; _) , (49) 
where normalization constant square is 

2 (\m± S \+n)l 

"< m±s nl(\m±s\\) 2 { ' 



and 



■J , 



/>OC 

/ e- p p\ m±s \ +k F 2 (-n,\m±s\ + l;p)dp (51) 
Jo 



The details of calculations of the more general integral can be found in [22 . 



00 



~- kz z u - 1 F 2 (-n,r,kz)dz= (52) 



e 







T(v)n\ I n ( n ~ — Pjyy — v — p — 1)(7 — v — p)... (7 — v + p) 



fc"7(7-l)...(7 + n-l) V 1 + -5 ((p+l)!) 2 7 (7 + l)-(7+p) 



Applying this formula to Eq. (|5Tj) one obtains 

n!(|m±s|!) 2 _ 1 
(\m±s\+n)\ ~ Cl m±s 



J n,± — n l „i , „m ~ 



Thus 



, n!(|m±s|+fc)!|m±s|! / ^ n\\m ± a\\ ngg(fc - p + q) \ 

n>± Qm±s\+n)\ \ ^ {{p+l)\y{n-p-l)\{\m± S \+p+\y: ) ' 1 ' 

fe= (|m + S | + fc)! / ^ ni \\ m + s\\Y{l^\k-p + q ) \ 

{\m-s\ + k)\ ( ^ n^m-sllllllVik-p + q) \ 

\m-s\\ y {{ P + iy.) 2 (n 2 -p-iy.(\m-s\+p + iy.j 1 ' 
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Particularly 



1\ = (|m + a\ + 1) ( 1 + 1 r") + (I™ - s l + !) i l + 1 —7 — t \ = 2n 



(57) 



I_ = \m + s\ — \m — s\ + 2(m — n 2 ) = 2n_ 



(58) 
(59) 
(60) 




l\ = (|m + s| + 3)(|m + s| + 2)(|m + s| + 1) + (|m - s| + 3)(|m - s| + 2)(|m - s\ + 1) + 
12 (m(|m + s\ + 3)(|m + s| + 2) + n 2 (|m - s| + 3)(|m - s| + 2)) + 
30 (m(ni - l)(|m - s| + 1) + n 2 (n 2 - l)(|m - s| + 1)) + 
20 (m(m - l)(m - 2) + n 2 (n 2 - l)(n 2 - 2)) 
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